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Abstract. The scope of the present paper is to determine how ion electrostatic 
wave perturbations in plasma flows are influenced by the presence of a 
kinematically complex velocity shear. For this purpose we consider a model 
based on the following set of physical equations: the equation of motion, the 
continuity equation and the Poisson equation for the electric potential governing 
the evolution of the system. After linearizing the equations, we solve them 
numerically. We find out that for a variety of specific values of parameters 
the system may exhibit quite interesting dynamic behaviour. In particular, 
we demonstrate that the system exhibits two different kinds of shear flow 
instabilities: (a) when the wave vectors evolve exponentially, the ion sound modes 
become unstable as well; while, (b) on the other hand, one can find areas in a 
parametric space where, when the wave vectors vary periodically, the physical 
system is subject to a strongly pronounced parametric instability. We also show 
the possibility of the generation of beat wave phenomena, characterized by a 
noteworthy quasi-periodic temporal behaviour. In the conclusion, we discuss 
the possible areas of applications and further directions of generalization of the 
presented work. 
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It is well known that plasma flows in different astrophysical, geophysical and 
laboratory situations are characterized by spatially inhomogeneous velocity fields 
(shear flows) and the presence of this velocity shear may significantly influence the 
modes of collective behaviour fostered in these flows. In most of the cases these flows 
are kinematically complex but even in a relatively simple case the overall dynamics, 
especially transition of the flow to turbulence, could be quite problematic [TJ [2] . One 
of the typical examples of astrophysical flows are helical plasma motions which occur 
in extragalactic jets [Sill] and in young stellar object jets |S]. It has also been argued 
that in the solar atmosphere at least some of macro-spicules are characterized by 
‘tornado-like’ kinematical flow geometry |6]. 

Obviously, the existence of such shear flows with kinematic complexity might 
strongly influence the plasma processes in a number of realistic astrophysical scenarios. 
Therefore, it is interesting to study the behaviour of plasmas being influenced by 
such kinematic complexity. In last two decades it has been realized that collective 
phenomena in shear flows are characterized by so-called non-modal processes, which 
in turn are related with non-normality of the involved mathematical operators |7]. 
In general, standard stability theory, i.e. normal mode analysis, does not describe 
completely the appearance of instabilities [ZIIHIIS]. An alternative approach is based 
on the method developed by Lord Kelvin m- In the framework of this method, the 
system of partial differential equations governing the dynamical evolution of modes 
of collective behavior is reduced to the inspection of a set of ordinary differential 
equations (ODEs) in time, i.e., to the solution of a relatively simple initial value 
problem. 

This approach can be effectively used in (magneto-)hydrodynamics, both for 
magnetized m M and unmagnetized plasma flows m- In particular, in |13j 
the authors considered the shear induced unstable modes applied to different, 
interesting astrophysical prototype structures - jets. The problem was examined 
for incompressible flows and it has been found that Alfven waves become subject 
to extremely strongly pronounced instabilities. It was shown that the shear flow 
instabilities may lead to the generation of large amplitude Alfven waves. A similar 
problem, for compressible magnetohydrodynamic flows, was studied by [I4j . where 
the authors have argued that the flow inseparably blends the slow and the fast 
magnetosonic and Alfven modes, leading to an efficient energy transfer from the 
background flow to the waves. 

In [15] the authors examined the electrostatic perturbations in an non-magnetized 
electron-ion plasma flow. Unlike the present paper, the authors considered the problem 
for a simple, one-dimensional, linearly sheared flow. It was shown that the ion-sound 
waves turn into plasma oscillations caused by a very efficient, shear-induced energy 
transfer mechanism between the mean flow and the waves. 

If the dissipation factors are efficient enough, the shear induced instability might 
lead to a substantial heating of the plasma flows nsun]. In particular, in m the 
authors considered acoustic waves and showed that the efficiency of the so-called self¬ 
heating by acoustic wave perturbations might be extremely efficient. The mentioned 
problem has been discussed in the context of non-magnetic chromospheric heating in 
solar-type stars. A very similar result has been obtained for magnetized flows m, 
where it was found that the rate of the self-heating mechanism might be high enough 
to be related with realistic heating scenarios in the solar atmosphere. 
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The scope of the present paper is to generalize the approach developed in m 
and to study the shear flow dynamics for ion electrostatic perturbations developing in 
more complex velocity configurations. In particular, as we have already emphasized, in 
previous studies the problem was considered for very simple, one-dimensional velocity 
shears. On the other hand, it is well known that in real astrophysical flows the 
kinematics might be quite complicated. Therefore, it is worthwhile to study more 
complex cases and see how the generalization, i.e. the increased degree of the flow 
complexity, alters the results qualitatively and/or quantitatively. In the present paper, 
we examine the problem physically, as a systematic plasma physics problem, without 
concrete astrophysical applications. Further applications to different astrophysical 
situations are in preparation and will be presented in separate publications in due 
time. 

The present paper is arranged in the following way. In the following section, we 
develop the theory of shear-induced instabilities for the ion perturbations. In the third 
section, we present and describe our results, while in the final section we summarize 
them and discuss them. 

2. Main Consideration 

In the present paper, we consider the non-magnetized, collisionless plasma flow with 
an electron temperature much higher than that of the ions, i.e. Tg 3> Ti. It is well- 
known that this kind of plasma sustains weakly damped low-frequency longitudinal 
electrostatic ion-sound waves with a constant ion acoustic speed 



where denotes the ion mass. 

We also assume that the quasi-neutrality condition holds, which means that the 
equilibrium electric field equals zero. On the other hand, the perturbations will 
inevitably lead to the generation of the perturbed electrostatic field, E = —Vc/, 
where (j) denotes the corresponding electric potential. Within the limits of the 
low-frequency approximation, the electron number density Ng is governed by the 
Boltzmann distribution 

Ng = noexp « no (^1 + • (2) 

The other basic equations governing the evolution of the system include the equation 
of mass conservation: 

DtN, + WV ■ V = 0, (3) 

the momentum conservation equation: 

AV = V(/., (4) 

rrii 

and the Poisson equation 

A(/ = 47re {Ng - Ni ), (5) 

where Ni denotes the ion number density, V denotes the flow velocity and Dt = 
+ (V • V) is the convective derivative. 
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In order to study the behavior of shear-induced instabilities, we linearize the 
system of equations around the equilibrium state: 

Ne,i = no + ne,i, (6) 

V = U + u, (7) 

where by Ue^i we denote the perturbed number densities of electrons and ions, 
respectively, and U and v are the unperturbed and perturbed flow velocity 
components. According to the linear approximation, it is assumed that all perturbed 
quantities are much smaller than the corresponding unperturbed quantities. After 
substituting Eqs. (15171) into Eqs. (EH they reduce to 

VtUi + NiV ■v = 0, (8) 

and 


Vtv + (v ■ V)U =-V(/>. 

rui 


(9) 


In accordance with the method developed in it is assumed that the 

unperturbed flow velocity U is spatially inhomogeneous. We then expand the velocity 
field in a Taylor series around the point A(xo,yoj zq), preserving only the linear terms: 

3 


U = U(A) 


1=1 


9U(A) 

dxi 


{Xj Xjo), 


with j = 1,2,3 and Xj = (x, y, z). 

One can straightforwardly show that the following ansatz 

F{x,y,z,t)=F{t)e^^-^\ 


with 


and 


3 

ijj-i ='^Kj{t)xj 

i=i 


3 „ 

V'z = ^Uj{A) I Kj{t)dt, 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


reduces the system of equations to a set of ordinary differential equations m- 
F(x,y, z,t) denotes the physical quantities n,Vx,Vy,Vz and F = {n,Vx,Vy,Vz} are 
the corresponding terms depending on time. By Kj{t) we denote the wave vector 
components, which are obeying the following differential equations m- 

atK + • K = 0, (14) 

where is a matrix transposed to the shear matrix S: 


S = 


/ Ux^x 

Ux,y 

Ux,z \ / All 

A21 

Ai3 \ 

1 

Uy,y 

Uy,^ = ^21 

A22 

A23 

V u.,x 

U.,y 

Uz,z ) \ Aoi 

A32 

A 33 / 


(15) 


where Ui^k = dUi/dxk- 

The plasma processes are described by Eqs. dSMl). Taking into account 
Eqs. (I11I12I13I) and omitting the symbol ”A”, one can rewrite them in the following 
dimensionless form: 

= k • u, 


(16) 
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4 ^) 

-h 

OllO^; 

+ ai 2 Uy -\- 

ax-iUz = 

-kx^, 


-h 

a 2 iUx 

a 22 Uy -\- 

023^2 = 

-ky(p, 

41) 

-h 

O-SlUx 

+ CL^lUy + 

03102 = 

-kzip, 


-f 

s'^ • k 

= 0, 




( 

On 

021 013 ' 

\ 


S = 


021 

022 O23 

h 



V 

031 

O32 O33 . 

/ 



d 


(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 


where d = inijn^^ u = vjCg 


= A,,/KMCs [hj = 1,2,3), ^ = -iTej)/e, 


oj = (jjp/Kz{0)Cs and tOp = ^^.-Kn^e^jmi corresponds to the ion-plasma frequency. 
Moreover, denotes the derivative of the function by a dimensionless time r = 
Kz{0)Cst, and Kz{0) corresponds to the initial value of K^- 

We intend to study the flow dynamics from the point of view of shear instabilities. 
For this purpose, in order to study the energy transfer related with the perturbations, 
we introduce their total energy: 


Etot = Y + Y’ 


(23) 


where the first and second terms are the kinetic energy and the compressional energy 
of the perturbations, respectively. 


3. Discussion 

In this section, we study numerically different interesting regimes related with the 
velocity-shear induced behavior of the perturbations. One can examine the problem 
by considering two different cases. Depending on the values of the shear matrix 
parameters, the wave vectors either evolve exponentially, or exhibit a stable character 
of evolution. 

This can be easily seen for the particular case: 013 = 023 = 033 = 0, 022 = —an, 
in which Eq. (I20p simplifies drastically so that kz becomes constant and the other two 
components of the wave vector obey the following set of equations: 

-r^k, + A,=0, (24) 

and 

ky^^^ -r^ky+Ay=0, (25) 

where F^ = ah -I- 012021 , A^ = -011031 — 021032 and Ay = 011 O 32 - 012 O 31 . 

It is clear, that when F^ > 0 the wave vector has an unstable temporal evolution, 
varying exponentially in time. This in turn, means that the physical system will 
inevitably undergo an unstable behavior. On the other hand, if F^ <0 the wave vector 
will depend on time periodically. But for specific values of the physical parameters, 
even in this case, the instability may arise in spite of the stable character of k(t). 

As a hrst example, we consider the situation when initially only the density 
is perturbed. In Fig. [1] we show the temporal evolution of the density and energy 
perturbations normalized by their initial values. Here, the following set of parameters 
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Figure 1. The temporal behavior of the normalized density dnir) = d(T)/d{0) 
and energy Enir) = Etotir)/Etot{0) perturbations. Here, the following set of 
parameters was chosen: a; = 1, an = —a 22 = 0.1, ai 2 = ai 3 = a 2 i = a 23 = 

0^31 = = ^^33 = 0, kxO = kyQ = 1, UxO = UyQ = UzO = 0, do = 0.1. 

was chosen: w = 1, an = -022 = 0.1, 012 = ais = 021 = 023 = 031 = 032 = 033 = 0, 
kxo = kyo = 1, Uxo = Uyo = U 20 = 0, do = 0.1, Therefore, T^ = 0.1 implying that 
the wave vector varies in time exponentially. In particular, as it is clear from the 
plots, for the considered time interval, r S [0 — 50], the amplitude of the density 
perturbations increases very rapidly from 1 to 80, which consequently leads to an 
increase of the energy of perturbations up to ^ 8 x 10^. Generally speaking, this 
means that the source of the amplification of the ion sound waves is the background 
flow energy, which in the framework of the present approach, behaves as an infinite 
energy reservoir. 

However, as the investigation of the shear flow dynamics shows, the system might 
undergo an extremely strong instability even for a negative value of T^. For studying 
this particular case, we examine the following set of parameters: an = —a 22 = 0, 
012 = —0.594, a 2 i = 0.2, while the rest of the parameter values are the same as 
in the previous case. It is evident that now F^ < 0 and, hence, the wave vector 
is characterized by periodic oscillations. In spite of this fact, the perturbations are 
strongly unstable. In Fig. [2] we present again the behavior of the density and energy 
perturbations. We see that the initially perturbed sound waves amplify very rapidly. 
In particular, the amplitudes of both the density and energy perturbations increase 
up to ^ 540 and ~ 2 x 10^, respectively. It is worthwhile to note that this instability 
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Figure 2. The temporal behavior of normalized density dnlr) = d{T)/d{0) 
and energy Enlj) = Etot(T)/Etot(0) perturbations. Here, the following set of 
parameters was chosen: oj = 1, ai 2 = —0.594, a 2 i = 0.2, an = aia = 022 = 
0,23 = “31 = “32 = “33 = 0, = kyQ = 1, u^q = UyQ = u^q = 0, dg = 0.1. Note 

that the range of ai 2 where the evolution of the modes is parametrically unstable, 
is very narrow, viz. [-0.607; -0.580]. 


disappears if one slightly changes the parameter values. As a matter of fact, for 
example, if we change (increase or decrease) 012 by no more than ~ 2.4%, the system 
becomes stable. One can straightforwardly check that the instability takes place 
only when ai 2 £ [—0.607;—0.580]. One can thus conclude that this is, as we have 
anticipated, an instability of parametric nature. 

In the hrst two examples above, we examined the evolution of some initially 
perturbed ion sound waves. These waves might be indirectly excited by velocity 
perturbations. This particular example is illustrated in Fig. [3l where the time 
evolution of the normalized total velocity, u„(t) = + Uyir^ + Uz{tY/ ug 

and density is presented. Here, the following set of parameters was chosen: a; = 1, 
Oil = 013 = O22 = 023 = O31 = 032 = 033, 012 = —0.594, O21 = 0.2, kxO = kyO = 1, 
UxQ = 0.1, Uyo = UzQ = 0, and do = 0. One can see from the upper graph that 
the normalized total velocity is strongly unstable, revealing the parametric instability. 
Even though initially only the velocity is perturbed, we see that in due course of time, 
density perturbations arise as well and, correspondingly, parametrically unstable ion 
sound waves are generated. 
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Figure 3. The temporal behavior of the normalized total velocity Unir) = 
u(t)/u{0) and density ^(t) perturbations. Here, the following set of parameters 
was chosen: ui = 1, an = ai 3 = a 22 = 023 = 03l = 032 = 033 , 012 = —0.594, 
021 = 0 . 2 , k^o = kyo = 1 , u^o = 0 . 1 , Uyo = UzO = 0 , do = 0 . 


It is well known that shear flows under favorable conditions exhibit so-called 
beat modes m- In Fig. 21 we show the behavior of both the normalized total 
velocity Unir) = u(t)/u(0) and energy Enir) perturbations. Here, the following set 
of parameters was chosen: w = 1, an = 013 = 022 = 0,23 = 032 = 033 , 012 = —0.594, 
021 = 0 . 2 , 031 = - 0 . 12 , = kyo = 1 , um = 0 . 1 , Uyo = u^o = 0 , and do = 0 . 

Figure [4] displays an example of such a process. From this figure we see that the 
ion sound wave ‘pulsates’, i.e. it is characterized by a quasi-periodic temporal “beat”. 
Such a structure is the result of the superposition of two frequencies that differ only 
by a very small amount. The fast oscillations are characterized by the ion acoustic 
frequency, that is modulated by a lower frequency defined by the shear parameters. 

4. Conclusions 

The goal of the present study was to consider ion-sound waves influenced by kinematic 
complexity and to study the role of the latter in the dynamical evolution of these waves. 
In particular, in the framework of the method developed by [10] we have considered 
the full set of equations, consisting of the momentum conservation equation, the 
mass conservation equation and the Poisson equation, and we linearized these around 
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Figure 4. The temporal behavior of normalized total velocity u„{t) = u{t)/u{0) 
and energy Enir) perturbations. Here, the following set of parameters was chosen: 
uj = 1, an = ai 3 = 022 = 023 = “32 = <133, ai 2 = —0.594, a 2 i = 0.2, 
a3i — 0.12, hxQ — tcyo — 1, — 0.1, UyQ — — 0, do — 6. 


an equilibrium state. We have shown that, depending on the choice of the set of 
parameters, the system might undergo a very efficient instability. In particular, it 
has been found that under favorable conditions, due to the velocity shear, the wave 
vector becomes unstable, making the ion-sound wave strongly unstable as well. On 
the other hand, we have also shown an example of the excitation of acoustic modes 
by means of only velocity perturbations. As yet another class of instability, we have 
found that even when the wave vectors behave periodically, for certain ranges of the 
parameter values, the system becomes strongly unstable. Moreover, it has been shown 
that the unstable character of the ion-sound waves dramatically changes to a steady 
behavior by slightly changing the parameter values. Another interesting feature of the 
examined system is that, under certain conditions, it exhibits an ”echo”-like behavior 
with quasi-periodic pulsations of the ion-electrostatic modes. 

In the near future, we plan to apply the developed method to realistic 
astrophysical flow models. In particular, it would be interesting to study the coupling 
of ion electrostatic waves with the velocity shear in stellar atmospheres, including 
young stellar objects. The idea that the ion-acoustic waves might influence the 
properties of stellar atmospheres has a long-standing story. For instance, in |18] 
the authors have considered the ion-sound waves in the solar atmosphere, studying 
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the transition of flat solitary waves into spherical modes. One of the intriguing 
problems concerning solar physics is the so-called chromospheric heating, that cannot 
be explained only by the convective component m- In the framework of acoustic 
waves, there have been proposed several mechanisms of chromospheric heating |201I21) . 
but observations with Transition Region And Coronal Explorer (TRACE) NASA space 
telescope have shown that there is about 90% deficit in the energy flux required to 
heat the chromosphere [22l [23l [24] . Therefore, it is interesting and quite reasonable to 
apply our model of kinematically driven ion-sound waves to the solar chromosphere 
and study the efficiency of heating. 

In the present paper, we studied the fluid in a non magnetized media, although in 
most of the astrophysical scenarios the magnetic field plays a crucial role. Therefore, 
it is of fundamental importance to generalize the present work by taking into account 
the appearance of cyclotron modes and their possible coupling via the agency of the 
velocity shear with ion-sound waves. 
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